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We study D- meson resonances in hot, dense quark matter within the NJL model and its Polyakov- 
loop extension. We show that the mass splitting between D^ and D~ mesons is moderate, not in 
excess of 100 MeV. When the decay channel into quasifree quarks opens (Mott effect) at densities 
above twice saturation density, the decay width reaches rapidly the value of 200 MeV which entails a 
spectral broadening sufficient to open J/'(/) dissociation processes. Contrary to results from hadronic 
mean-field theories, the chiral quark model does not support the scenario of a dropping _D-meson 
masses so that scenarios for J /^l> dissociation by quark rearrangement built on the lowering of the 
threshold for this process in a hot and dense medium have to be reconsidered and should account 
for the spectral broadening. 

PACS numbers: 12.39.-x, 12. 38. Aw, ll.lO.St, 12.38.Lg 



I. INTRODUCTION 

The modification of the D meson properties (masses 
and widths) in hot, dense matter has consequences for 
scenarios of J/V' suppression, e.g., by the processes of 
the type 

J/iJj + TT^p^ D* +D , J/iP + N ^ Ac + D , (1) 

which couple hidden charm to open charm states and 
thus lead to a dissociation of charm in the medium [l| , see 
Refs. |3t3 for early controversial estimates of the cross 
sections of such processes. The reverse process @-13 of 
charmonium regeneration by open charm recombination 
should play a dominant role for J/ip production at LHC 
[9|, llO| where charm is abundant in the medium. 

Since either dropping masses [ll|, [ij] or increasing 
widths [Ij, [ij] of the D mesons in a hot and dense 
medium can lead to a lowering of the reaction thresh- 
old and thus to a strong increase of the rate for the pro- 
cesses ([1]) , both effects may contribute to an explanation 
of the anomalous J/ip suppression found in the NA50 ex- 
periment [iSHlTI I and subsequently confirmed by NA60 
[ll[ii] and PHENIX ^. For a recent review, see [2l[. 

In contrast to results from a relativistic mean-field 
model of D mesons in nuclear matter which predicts a 
strong downwards shift of £)+ meson masses due to the 
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renormalization with a scalar mean field [22| , the consid- 
eration of the quark substructure of these mesons leads 
to qualitatively different behavior. As we will show in 
this work on the basis of a chiral quark model of the 
Nambu-Jona-Lasinio (NJL) type and its Polyakov-loop 
extension (PNJL), the Pauli blocking effect in the Bethe- 
Salpeter equation for the L'-mesons largely compensates 
the dropping masses of their quark constituents. As a 
result, D-meson masses do not drop but stay almost 
constant or rather increase with increasing density (and 
temperature) of the matter. Their decay width, how- 
ever, increases rapidly and reaches values which allow 
for a subthreshold quark rearrangement dissociation of 
J/ip. Therefore, scenarios for J ftp suppression built on 
the quark rearrangement reaction ([T|) have to be recon- 
sidered. As has been demonstrated in ^,i2J], a sufficient 
width of D mesonic correlations in the quark plasma is 
essential for understanding charm thermalisation and dif- 
fusion in RHIC experiments, see [23 for a review. 

It is interesting to note that a recent self-consistent 
coupled channel a ppro ach for D mesons in hot, dense 
nuclear matter [26|, |27| supports the picture of a spec- 
tral broadening with a negligible mass shift up to tem- 
peratures T = 150 MeV and densities n = 2 uq with 
no = 0.16 fm~'^ being the nucleon density of nuclear 
matter at saturation. Therefore, it seems likely that a 
quark hadron duality similar to that discussed for low- 
mass dilepton production [2^ can be observed also in the 
D meson channel at the deconfinement transition. The 
spectral broadening of D mesons rather than their mass 
shift has been suggested for an explanation of anoma- 
lous J/tp suppression in [l3lll4l|. This question, however, 
awaits a thorough investigation. 



In the present note, we investigate the extension of 
previous exploratory calculations of D mesonic correla- 
tions in quark matter, based on the NJL model [lj,[2al to 
the domain of finite baryon densities which will become 
accessible in the CBM experiment at FAIR. Of partic- 
ular interest will be the question whether the suggested 
mass splitting of D meson states [30| will be observable or 
rather washed out by spectral broadening. Furthermore, 
a strong isospin dependence of the D meson broadening 
could result in observable signatures, possibly relevant 
for quark-gluon plasma diagnostics. 

Chiral dynamics has been applied successfully not only 
in the light quark sector but also especially for the in- 
vestigation of heavy-light pseudoscalar meson properties. 
This has been most impressively demonstrated within the 
Dyson-Schwinger equation approach in Ref. [3l| which 
reproduced a complete set of heavy meson observables. 
A particular feature of these systems is the heavy quark 
symmetry which allows to separate the physics of the 
heavy and the light quark components from each other 
and absorb the chiral dynamics of the latter into the uni- 
versal Isgur-Wise function j^. For more details, see 
the reviews on heavy quark effective theory [33l - l35| . It 
is interesting to note that already the chiral dynamics 
encoded in the rather schematic NJL model reproduces 
features like heavy-quark symmetry and Isgur-Wise func- 
tion [36| when extended to the heavy quark sector. These 
properties remain unaffected when confining properties 
are mimicked in the heavy-quark extended NJL model 
by an infrared cutoff procedure |37| or by a confining 
interquark potential in a rclativistic potential model of 
heavy mesons [SSj . 

D meson properties as reported in our study come from 
a simultaneous solution of two types of nonperturbative 
equations: the gap equations for the quark masses and 
the Bethe-Salpeter equations for the meson masses. The 
results for the dependence of pseudoscalar meson masses 
Mp on the current quark masses mq of the heavier quark 
in the meson nicely reproduce the transition from a Gcll- 
Mann-Oakes-Renner like behaviour for pion and kaon 
{Mp ^ J'niq) to the additive quark model like behaviour 
{Mp ^ mq) for D and B mesons. These results are 
then readily generalized from the vacuum to finite tem- 
perature and chemical potential within the Matsubara 
formalism. 



II. MODEL AND FORMALISM 

We employ a four - flavor model with NJL - type inter- 
action kernel as a straightforward generalization of recent 
work on the SU/(3) scalar and pseudoscalar meson spec- 
trum |39l - l41j developed on the basis of Ref. ^^ and its 
generalization by coupling to the Polyakov loop [4^, |4J1 , 

15 

£ = q{irD^ + m)q + GsY.[{q\''qf + {qil5X''qf 

a=0 

- U {'^[AmA]\T) . (2) 



Here q denotes the quark field with four flavors, Nf = 4, 
/ = u,d,s,c, and three colors, Nc = 3; A° are the fla- 
vor SU/(4) GeU - Mann matrices (a = 0, 1, 2, . . . , 15), 
Gs is a coupling constant. The global symmetry of the 
Lagrangian ([2]) is explicitly broken by the current quark 
masses m = diag(?7i,„, TTi^, ttIs, "nic)- The covariant deriva- 
tive is defined as D^^ = 9'^ - iA^' , with A^' = (5^^" 
(Polyakov gauge); in Euclidean notation A^ — —iA^. 
The strong coupling constant gs is absorbed in the defi- 
nition of A^'{x) = gsAi^{x)^, where Ai^ is the (SUc(3)) 
gauge field and Aq are the (color) Gell-Mann matrices. 

The Polyakov loop field $ appearing in the potential 
term of ([2]) is related to the gauge field through the gauge 
covariant average of the Polyakov line [45| 



$(f) = ((Z(x))) = ^Tr,((L(x))), 



where 



L{x) = Pexp 



drAji^ {x, t) 



(3) 



(4) 



Concerning the effective potential for the (complex) $ 
field, we adopt the form and parametrization proposed 
in Ref. 0. 

This effective chiral field theory has the same chiral 
symmetry of QCD, which is also shared by the quark 
interaction terms. The (P)NJL model is a primer for 
describing the dynamical breakdown of this symmetry 
in the vacuum and its partial restoration at high tem- 
peratures and chemical potentials. At the same time it 
provides a field-theoretic description of pseudoscalar me- 
son properties which is in accordance with the low energy 
theorems (such as the Goldstone theorem) of QGD. 

In the vacuum the PNJL model with the Lagrangian 
([2]) goes over to the NJL one and the pseudoscalar meson 
properties arc described in the standard way by analyzing 
the polarization operators 

n.,(P) = ^Nc J -0^tTD [S.{p){^^,)S,{p + P)(z75)] ,(5) 

where trn is the trace over Dirac matrices, Si{p) is the 
quark Green function with the dynamical quark mass Mi . 
The polarization operators can be presented in terms 
of two integrals which for mesons at rest in the medium 
are given by 



-[P^-{M,-M,f]l^{Po)}, (6) 



where 



II = iN, 



d^p 



Nr. 



A „2 



P cfp 



Il'iPo) 



47r2 



(27r)4 pI - Ej 47r2 J, E, 
d^p 1 



(7) 



(2^)4 {pl-EmP^, + P,Y-E^) 






E, + E, 



P2 



(^^ 



■E,Y 



(8) 



where Eij = J-p^ + Mf- is the quark energy. 

As the Lagrangian ^ defines a nonrenormahzable field 
theory, we introduce the 3 - momentum eutoff with the 
parameter A to regularize the integrals. When Pq > Mi + 
Mj it is neeessary to take into aceount the imaginary part 
of the second integral. It may be found, with help of the 
Pq — ie, that 



ie -prescription Pq 



Il\P^) 



Nc, 
4^' 



A 



)2dp 



E, + E, 



Nr. 



^i^j Pq - i^i 
P* 



E,r 



167r {E* + E 



(9) 



(M, 



M,)2)/2Po 



where p* = ^{P-j - (M, - M,y){P,' - 
is the momentum and E*j ~ \ {v*Y 
sponding energy. 

The quark mass Mi we find from the gap equation 



Mf, the eorre- 

''iJ 



M, 



16M,GsPi. 



(10) 



The meson mass spectrum we obtain from the condition 



1 - 2Gsn*^(Po = Mp, P = 0) = 



(11) 



The pseudoscalar meson-quark-antiquark coupling con- 
stants are defined as 



-2 



1 



d 



2Mp dPn 



[n*-'-(Fo)] 



=Mp 



(12) 



Note that when Pq > J^h + Mj, then Eq. ([TT]) has to 
be solved in their complex form in order to determine 
the mass of the resonance Mp and the respective decay 
width Tp. Thus, we assume that this equation can be 
written as a system of two coupled equations 



Mf, - ^Tl 



(A/, - Af,f 



{8Gs)-'-in+Ii) 
lUPo ^ Mp + ie)\^ 
- [MpTp] = 

{8Gs)-'-{n + Ii) 
MPo = Mp + ie)\^ 



Rel2{Po = Mp + ie), (13) 



liJihiPo^Mp + ie), (14) 



which have solutions of the form 



Po = Mp 



1 
^^Tp. 



(15) 



As shown in [30], the model ^ suecesfuUy describes 
meson properties in the vacuum at T = /i = 0. We use 
here the paramctrization given in table III of Ref. [43| 
for the case of the NJL model. The value of the coupling 
constant is GgA^ = 2.32, the three-momentum cutoff 
is at A = 602.3 MeV, see also [4^ for an online tool 
fo the three-flavor NJL model. The parametrization of 
the current-quark masses for the heavy flavors (c, b) is 
performed with the above formulae for the masses of the 
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FIG. 1: (Color online) The systematics of heavy-light pseu- 
doscalar meson masses Mp in dependence on the mass m/ 
of the heavy flavor in the bound state with a light antiquark 
within the present model (red solid line) is compared to the 
masses of heavy-light pseudoscalar mesons (tt, K, D, B) and 
the limits for the corresponding current quark masses (sym- 
bols with error bars) according to the particle data group [401. 



corresponding heavy-light pseudoscalar mesons {D,B). 
The results are summarized in table HI The dependence 
of the heavy-light meson mass on the current mass of the 
heavier quark is shown by the solid line in Fig. [1] and 
provides satisfactory agreement with the particle data 
group listings [49| . 



flavor 


p 


Mp fp 


771/ Mf 






[MeV] [MeV] 


[MeV] [MeV] 


u,d 


TV 


135.0 92.4 


5.5 368.0 


s 


K 


497.7 95.4 


140.7 587.4 


c 


D 


1869.3 79.6 


1279.9 1828.8 


b 


B 


5279.4 


4634.8 



TABLE I: Results for pseudoscalar meson properties in the 
light, strange, charm and bottom sectors for the correspond- 
ing current quark masses and dynamically generated quark 
masses of the model in the vacuum at T = /x = 0. 



The generalization of the model ([2|) to the case of 
nonzero temperature and density (details in [39| ) is 
done within the imaginary time formalism by introduc- 
ing the Matsubara frequencies w„ = (2n -I- l)7rr, n = 
0, ±1, ±2, . . ., so that Po — > ii^n + M with the chemical 
potential /z and the temperature T . Instead of the in- 
tegration over poi we have now to perform a sum over 
Matsubara frequencies. In the result we obtain the inte- 



grals 



1 47r2 



E, 



(16) 



-Nr 



1 






2^, (£;, + Po)2 - ^' 



1 



2EUE, - Po)^ - E^ ' 
1 1 + 

Tl 

2E, {E, ~ PoY ~ Ef ' 

1 1 

'2E,{E,+Poy-Ef "^■ 



(17) 



where n^ = /$(±i?i) are the generahzed fermion distri- 
bution functions [4J] for quarks of flavor i with positive 
(negative) energies in the presence of the Polyakov loop 
values $ and $ 



ME) = 



1 + 3(1, + $g-/3(£;-/x))g-/3(£;-Ai) _^ g-3/3(B-^0 ' 



(18) 



which go over to the ordinary Fermi functions in the case 
of the NJL model, where $ = $ = 1 



fii±E) = 



1 



1 + QP{±E^t,} 



(19) 



Note that we put the chemical potential for charm quarks 
to zero in the calculations discussed below. A small finite 
value might be necessary to ensure exactly vanishing net 
charm at high baryon density when, as we demonstrate 
below, the symmetry between masses of D mesons with 
charm and those with anticharm is broken by medium 
effects. 



III. RESULTS 

We have performed selfconsistent solutions of the gap 
equations for the dynamically generated light quark 
masses at finite temperatures and chemical potentials for 
the PNJL model and its NJL model limit within the stan- 
dard setting as summarized above. In Fig. [2] we display 
the results for the restoration of the approximate chiral 
symmetry in the u, d quark sector along trajectories in 
the QCD phase diagram with a constant ratio r = T/ fi, 
where r = 0, 1/3, 1/2, 1. The (pseudo)critical tem- 
peratures for the NJL model and for the PNJL model 
are shown in the phase diagram in Fig. [3] together with 
the trajectories along which we investigate the D meson 
properties. 

When lowering the ratio r — > 0, the phase transition 
turns from crossover to first order. The chiral restoration 
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FIG. 2: (Color online) Dynamically generated light quark 
mass Mu(d) as a function of the quark chemical potential along 
trajectories in the quark matter phase diagram for T = 
(black solid line), T = /i/3 (green dotted line), T = fj,/2 (blue 
dash-dotted line) and T = jj, (red dashed line). Results of the 
NJL model (thin lines) are compared to those of the PNJL 
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FIG. 3: (Color online) (Pseudo)critical temperatures versus 
quark chemical potential for the NJL model and for the PNJL 
model. The first order transition lines are shown as solid lines, 
the crossover is indicated by dashed lines. Critical end points 
are denoted as full dots. The straight lines stand for the 
trajectories T — rjj, along which we investigate the D meson 
properties: r — 1, 1/2, 1/3. 



is a result of the phase space occupation (Pauli blocking) 
which effectively reduces the interaction strength in the 
gap equation. 

Naively, one would expect that heavy-light mesons 
such as D mesons, should suffer a mass reduction when 
embedded in a hot and dense medium, as a result of the 
strong reduction of the light quark constituent mass (the 
charm quark mass is approximately unaffected). The so- 
lution of the in-medium D~ meson Bethc-Salpcter equa- 
tion (BSE), however, shows a different result, displayed 



in Figs. I4I5I for the NJL case. This can be understood 
from inspecting the kernel of the BSE: The mass shift of 
the light quarks, which affects the energy denominators 
and would lead to a lowering of the meson masses, when 
evaluated in free space, gets compensated or even over- 
compensated by the Pauli blocking factor in the numer- 
ator. It is clear that D mesons containing a light quark, 
as the dominant quark species in a dense medium, feel a 
stronger Pauli blockig than those containing light anti- 
quarks. Therefore, the D~ mass is shifted upwards while 
the D^ mass stays approximately constant. 
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Polyakov loop. The results are rather similar, except for 
the smaller width of the transition region in the PNJL 
case and the different position of the critical point in the 
phase diagram. The higher temperature of the critical 
endpoint entails that along the trajectory for r — 1/3 
the medium undergoes a first order transition and the D 
meson properties change discontinuously. 



FIG. 4: (Color online) Density dependence of D meson masses 
(upper panel) and widths (lower panel) along the trajectory 
T = fi in the QCD phase diagram for the NJL model case. 
While the D meson bound state masses are almost constant or 
moderately rise with increasing density the continuum thresh- 
old (M„{rf) + Mc) gets lowered dramatically due to the chiral 
symmetry restoration transition. At about twice nuclear den- 
sity, the decay channel into their quark-antiquark constituents 
opens and the bound states become resonant scattering states 
in the quark-antiquark continuum (Mott effect). 

There is an important consequence of this fact that 
the heavy-light continuum threshold is lowered while the 
bound state masses are not: at a critical density the 
bound states merge the continuum (Mott effect) and 
become unstable against decay into their quark con- 
stituents, as signalled by the nonvanishing decay widths 
shown in the lower panels of Figs. JH |5l The stronger 
Pauli blocking for the D~ mesons leads not only to a 
lower Mott density, but also to a larger decay width as 
compared to the Z?+. From these Figures one can also 
read off the critical densities where the phenomenon of D 
meson Mott effect could be expected for the NJL model. 

In Figs. I6I7I we present results as a function of the 
quark chemical potential in isospin symmetric matter for 
the NJL model and its extension with coupling to the 
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FIG. 6: (Color online) D meson masses and widths along the 
trajectories T — n, T = fi/2 and T = n/S in the QCD phase 
diagram for the NJL model. 



The modification of D meson properties in hot and 
dense nuclear matter as reported in the present work is 
essentially different from the one suggested in ReL 



II 



where a lowering of the D meson mass had been con- 
jectured with consequences for charmonium dissociation 
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FIG. 7: (Color online) Same as Fig. [6] for the PNJL model. 



in heavy-ion collision experiments. In this reference the 
Pauli blocking effect was neglected. It is interesting to 
note that the Pauli blocking effect occurs not only on the 
quark level but also on the hadronic level of description, 
when coupled channel equations for D mesons in nuclear 
matter are solved selfconsistently as, e.g., in Ref. [26| . 
Also in this approach the D meson mass remains almost 
unshiftcd while a considerable spectral broadening is ob- 
tained under similar conditions of density and tempera- 
ture as considered in the present work. 



quarks (£)+ — dc, D° = uc) suffer no Pauli shift since 
there are no antiquarks in the medium. Their mass at 
the Mott transition density is approximately the same as 
in vacuum. 

It is interesting to compare the decay widths of the D 
mesons into their quark constituents. Due to the repul- 
sive Pauli shift, the D~ and Z)° are clearly above the 
threshold already at the first order chiral phase transi- 
tion and have a non-negligible decay width, whereas the 
D^ and D^ mesons are still very good resonances with 
negligible decay width at the transition. 

The effect of coupling the chiral quark model to the 
Polyakov loop in the PNJL model is an effective sup- 
pression of the quark distribution functions in the BSE 
(fT51) . (fHl) as long as $ ^ I which leads to a narrowing 
[5 of of the chiral transition region where medium effects 
on the D meson properties occur in the present model. In 
this region the dissociation of D mesons (Mott effect) oc- 
curs and is signalled by an increase in the spectral width 
of these states. 

Summarizing the results of this model, we find no sup- 
port for dropping D meson masses in the vicinity of the 
chiral/deconfinement phase transition in hot and dense 
matter but rather a strong spectral broadening. There- 
fore, scenarios of J/ip suppression in dense matter via dis- 
sociation processes like ([T]), which are built on increasing 
widths [ij, ^^ for the D mesons should be conceptually 
preferable over those built on dropping masses. 



IV. CONCLUSIONS 

In the isospin-symmetric quark matter case holds 
M£,+ = Mijo and M^- = Mqo . 

The D mesons containing light quarks {D^ = cd, 
Z)° = cu) suffer a positive mass shift due to the effec- 
tive reduction of the coupling by Pauli blocking, since 
the phase space is occupied by light quarks abundant 
in the medium. The D mesons composed of light anti- 
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